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Abstract 

Fast heuristically weighted, or pseudo-C^, estimators are a frequently used method 
for estimating power spectra in CMB surveys with large numbers of pixels. Recently, 
Challinor & Chon showed that the E-B mixing in these estimators can become a 
dominant contaminant at low noise levels, ultimately limiting the gravity wave signal 
which can be detected on a finite patch of sky. We define a modified version of the 
estimators which eliminates E-B mixing and is near-optimal at all noise levels. 
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1 Introduction 



Measurements of CMB polarization, reported by four experiments to date 
fjjtllll HI Hob . are rapidly improving in sensitivity. Future ground-based surveys 
will be optimized for sensitivity to the B-mode power spectrum, which offers 
a unique window into the physics of the early universe via the gravity wave 
sign al (flit 0) , and can help break parameter degeneracies via the lensing signal 



On a practical level, two general methods for estimating power spectra from 
noisy maps have been used in experiments to date: maximum likelihood (0), 
and pseudo-C^ estimators (flU |(| HI). The strengths and weaknesses of these 
methods are complementary. Maximum likelihood power spectrum estimata- 
tion is optimal (in the sense that the Cramer- Rao inequality is saturated) but 
computationally expensive, requiring 0(N^ ix ) time and 0(iVp ix ) memory, and 
rapidly becomes infeasible for large maps. In contrast, pseudo-CV estimators 
are very fast, suboptimal in principle, but frequently near-optimal in practice. 
For this reason, they are currently the method of choice for surveys with large 
pixel counts. 
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In the context of CMB polarization, pseudo-CV estimators have an additional 
drawback, as stndied by Challinor & Chon fl" the estimated BB power spec- 
trum acquires a nonzero contribution from E-modes in the map. The pseudo-C^ 
construction includes a debiasing step (Eq. (5) below) which removes the E-B 
mixing in the mean; however, the variance of the BB estimators still depends 
on the level of EE power. This is analagous to treating noise by subtracting 
the noise bias from each estimator; even though the bias from noise is removed 
in the mean, it still makes a contribution to the estimator variance. This con- 
tribution can dominate if the instrumental noise is sufficiently small; e.g. in 
(0) it is shown that for surveys with f s ^ y ~ 0.01, it limits the gravity wave 
signal which can be detected using pseudo-CV estimators to T/S ~ 0.05. 

We will describe a modification to the pseudo-C^ construction, defining pure 
pseudo-Ci estimators which separate E and B in the strongest possible sense. 
On a finite patch of sky, the BB estimator completely filters out E-modes, 
acquiring contributions only from the B-mode signal and noise. We will show 
that these estimators significantly improve power spectrum errors for noise 
levels < 20 /iK-arcmin, and impose no limit on the gravity wave signal which 
can be detected. 



2 The pseudo-CV construction 



In this section, we briefly review the pseudo-CV construction before presenting 
our modification in §3. For more details, see (|5j). 

The basic idea is easy to explain. One constructs power spectrum estimators 
on a finite patch of sky by weighting the observed polarization map Tl a b(x) by 
a heuristically chosen weight function W(x), then computing a^ m 's and power 
spectra as if the weighted polarization were an all-sky field: 



a E im A =\ d 2 x2W(x)U ab (x)Y ( ^ )ab (x) (1) 
~aV=\ d 2 x2W(x)U ab (x)Y^ )ab (x) (2) 

^ffffdef I S^„E*„E . p^BB def 1 „B* n B (o\ 

~ 2£ + 1 im a im W _ 2£ + 1 171 m 

where Y^, denote E-mode and B-mode spherical harmonics Klh. We 
have introduced the notation Ci for the power spectrum of the weighted po- 
larization field, or "pseudo power spectrum". Because multiplication by the 
weight function in real space mixes multipoles in harmonic space, the pseudo 
power spectrum Ci acquires contributions from multipoles £' ^ £. In the mean, 
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this contribution can be written in terms of transfer matrices K^,: 




Kz K + 




N e EE 



(4) 



where the vectors Ng represent noise bias. Note that the matrices K^,, K^, 
represent multipole mixing from the weight function with and without accom- 
panying mixing of E and B. 

The final ingredient in the pseudo-CV construction is a debiasing step: 

1 (, 




K~ K + 



C EE - N EE 



\ 



C BB - No, 



(5) 



This removes the bias imposed by the weight function and noise; each Cg is an 
unbiased power spectrum estimator. (As a technical point, we mention that 
the matrix inversion in Eq. (5) can only be performed if / s k y is large; otherwise 
one has to bin multipoles into bandpowers.) 



3 E-B mixing 



An elegant framework for studying E-B mixing has been presented in (|3j), 
which was the original inspiration for this work. On a finite patch of sky, 
the polarization field can be decomposed into three types of modes: pure E- 
modes, which receive no contribution from the B-mode signal, pure B-modes, 
which likewise receive no contribution from the E-mode signal, and ambiguous 
modes, which receive contributions from both. 

In terms of this decomposition, the underlying reason why pseudo-CV estima- 
tors mix E and B modes on a finite patch of sky can be understood as follows. 
Returning to the definition in Eq. (2), the pseudo multipole af m is given by 
the overlap integral with a mode (in brackets below) which is a mixture of all 
three types: 



af m = I d 2 x2U ab (x) 



W{x)Y£(x) (6) 



Therefore, the E-mode signal makes a nonzero contribution to af m , via the pure 
E-mode and ambiguous mode components of the mode in brackets. In this way, 
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the E-mode signal finds its way into the BB power spectrum estimator Cf B , 
where it contributes extra variance. 

Ideally, one would like to have a BB power spectrum estimator which does 
not receive contaminating contributions from the larger E-mode signal. This 
suggests modifying the definition in Eq. (6) so that the mode in brackets is a 
pure B-mode. 

According to (0), pure B-modes on a finite patch of sky are of the form 

Qe ac V c V fe + le fec V c V a ) $(x) (7) 



where <&(x) is a potential which must satisfy the boundary conditions $ = 
V a $ = 0. To make contact with the pseudo multipole in Eq. (6), we write the 
B-mode spherical harmonic as the result of applying the same differential 
operator to the scalar spherical harmonic Y^ m : 



l lm 



N e I d 2 x2U ab (x) 



W(x) -e ac V c V fc + -e fec V c V a Y t 



lm 



X 



where N £ = 1/^-1)^+1)^ + 2). Comparing Eqs. (7) and (8), it is nat- 
ural to modify the definition of af m by bringing the weight function inside 
the differential operator. We take this as the definition of the pure pseudo 
multipole: 



a 



B,pure def 
lm 



d 2 x2Tl ab \ 



x 



\e ac V c V b + ^6 bc V c V a ) W(x)Y* n 



x) 



(9) 



For the mode in brackets to be a pure B-mode, the weight function must 
satisfy the boundary conditions W = V a W = 0. We will return to this in §4. 

We define pure pseudo-Ci estimators by replacing df m af r £ ure and leaving 
subsequent steps in the pseudo-C^ construction from §2 unchanged. When this 
modification is made, the E — > B term in the transfer matrix is always zero: 



KZ) 



+ ,pure 



(10) 



This corresponds to the statement that pure pseudo-CV estimators do not mix 
E — > B in the mean. The estimators also separate E and B in a much stronger 
sense: in a single realization, the estimated BB power receives zero contribu- 
tion from E-modes in the map. This is because the observed polarization Il ab 
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always overlaps a pure B-mode (in brackets in Eq. (9)) which by construc- 
tion completely filters out E-modes, even on a finite patch of sky. Because of 
this strong separation, the pure pseudo-C^ estimator is completely "blind" to 
E-modes and EE power does not act as an source of extra variance. 

An essential ingredient in the pseudo-C^ construction is fast (C(^m ax )) algo- 
rithms for evaluating the estimators and precomputing the transfer matrix. 
Variants of these algorithms can also be given for pure pseudo-C^ estimators 
fflih . but we omit the details here. 



4 Statistical weight 

An important practical issue, for both pure and ordinary pseudo-CV estimators, 
is choosing the pixel weight function W(x). In general, this must be done 
empirically (perhaps guided by heuristics) using Monte Carlo simulations to 
optimize the estimator variance. The optimal weight function will depend on 
£; at low i or high signal-to-noise, uniform weighting (which minimizes sample 
variance) is near-optimal, whereas at high i or low signal-to-noise, inverse 
noise weighting (which minimizes noise variance) is near-optimal (0). 

For pure pseudo-C^ estimators, an extra complication arises when choosing the 
weight function. The statistical weight of a pixel is given by a combination of 
the weight function W(x) and its first two derivatives in the pixel. This can 
be seen from Eq. (9), where the estimator is defined by overlapping the map 
Tl ab (x) with a pure B-mode (in brackets) which contains terms proportional to 
W, its first derivative, and second derivative. In contrast, for ordinary pseudo- 
Ci estimators, the map overlaps a mode which is proportional to W(x), and 
the statistical weight of a pixel is simply given by the weight function. 

As a concrete example (taken from Ref. (flih). the weight function W(x) shown 
in Fig. 1, left panel, is optimized for pure pseudo-CV power spectrum estima- 
tion at i = 30 on a circular patch of sky with radius 13° and uniform white 
noise. If this weight function were used in ordinary pseudo-CV estimation, the 
statistical weight would be concentrated near the center of the survey and 
the resulting estimator would be suboptimal (since the noise distribution is 
assumed uniform). For pure pseudo-C^ estimators, one can understand the 
distribution of statistical weight by plotting the pure B-mode 



which overlaps the map U ab (x) in Eq. (9). This is shown, for (£,m) = (30,0), 
in the right panel of Fig. 1. It is seen that the statistical weight is distributed 




(11) 
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Fig. 1. Weight function W(x), satisfying the boundary conditions W = V a W = 
required for pure pseudo-Q power spectrum estimation (left panel) and the resulting 
pure B-mode defined by Eq. (11) (right panel). 

roughly evenly throughout the survey. (For the m = mode shown here, 
the statistical weight is somewhat biased toward the center of the region, but 
would be uniformly distributed after averaging over m as in Eq. (3).) 

As remarked in §3, a general feature of pure pseudo-CV estimation is that 
the weight function must be apodized so that the boundary conditions W = 
VaW = are satisfied. Suppose that one tried to "cheat" by using a weight 
function which is constant throughout the survey except for a tiny ribbon 
around the boundary, where it rapidly goes to zero in order to satisfy the 
boundary conditions. Then the resulting estimator would have most of its 
statistical weight in the ribbon, where the derivatives become large, and would 
therefore be very suboptimal. In fact, one can think of the boundary conditions 
W = V a W = as a consequence of the statistical weight being given by the 
first two derivatives. If the boundary conditions were not satisfied, then the 
estimator still makes sense mathematically but the statistical weight would 
include delta functions on the boundary, and does not make sense for noisy 
data. 

Let us emphasize that, although pure pseudo-C^ estimators require the weight 
function W(x) to be apodized near the boundary, the distribution of statis- 
tical weight need not be apodized. Any distribution of statistical weight can 
be obtained by choosing W(x) such that the boundary conditions are satis- 
fied, and the amplitude of the pure B-mode in Eq. (11) matches the desired 
distribution, For example, the weight function considered in this section has 
been engineered to give uniform statistical weight (to match the assumption 
of uniform noise). This example is part of a systematic framework for opti- 
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Fig. 2. BB power spectrum errors for pure pseudo-C^ estimators (blue/left), ordi- 
nary pseudo-C^ estimators (green/middle), and the Cramer-Rao bound (red/right). 
In the left panel, we have shown a model with T/S = 0.2 and noise level 20 
/iK-arcmin; in the right panel, a model with T/S = 0.05 and noise 10 /iK-arcmin. 
The dashed lines show the lensing component of the BB power spectrum; detecting 
the gravity wave signal requires measuring power in excess of this level. 

mizing weight functions for both pure and ordinary pseudo-C^, which will be 
presented separately (flit). 



5 Examples 



We now consider some results for a mock survey on a spherical cap with radius 
13°, with uniform white noise and a Gaussian beam with #fwhm = 25 arcmin. 
In Fig. 2, we compare the power spectrum errors obtained using pure pseudo- 
Cg estimators, unmodified pseudo-CV estimators, and the Cramer-Rao bound 
on the estimator variance. 

It is seen that the pure pseudo-CV construction significantly improves power 
spectrum errors, relative to unmodified pseudo-CV, for noise levels < 20 /xK- 
arcmin. For these low noise levels, the extra variance in unmodified pseudo-C^ 
estimators can limit the gravity wave signal which can be detected. For exam- 
ple, in the right panel (10 /iK-arcmin), a gravity wave signal with T/S = 0.05 
can be detected using pure pseudo-Q estimators, but not using the unmodified 
versions. 

In Fig. 3, we have shown the values of T / S which can be detected at la using 
each estimator, for a variety of noise levels. We find that the smallest gravity 
wave signal which can be detected using unmodified pseudo-Q estimators is 
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Fig. 3. Minimum T/S detectable at la, for a spherical cap shaped survey with 
radius 13°, uniform white noise, and 25 arcmin beam. 

T/S = 0.042, in agreement with Challinor & Chon (4), but that the pure 
versions are ~ 80% optimal all the way to the floor at T/S ~ 10~ 3 , which 
represents the limit from the lensing component of the BB power spectrum 
(treated as a Gaussian contaminant) for this survey region. 



6 Discussion 

For CMB surveys with many pixels, pseudo-C^ power spectrum estimation 
is currently the "industry standard" , but has a shortcoming for future high- 
sensitivity polarization experiments: the EE signal acts as an extra source of 
noise on the BB power spectrum, becoming a significant contaminant at noise 
levels < 20 /zK-arcmin. We have defined pure pseudo-C^ estimators, which 
eliminate this shortcoming while preserving the fast (C(^inax)) evaiua tion of 
the estimator and transfer matrix. These estimators completely filter out the 
E-mode component of the map; the estimated BB power receives contributions 
only from the B-mode part of the signal and from noise. 

A property of the pure pseudo-C^ estimator is that the statistical weight is 
given by a combination of the weight function W(x) and its first two deriva- 
tives. This imposes no restriction on the way the statistical weight is dis- 
tributed (one can always choose W(x) to give any desired statistical weight 
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while satisfying the boundary conditions) but makes optimizing the weight 
function for a given noise distribution less intuitive. The optimization prob- 
lem will be studied systematically in (18), in the context of both pure and 
ordinary pseudo-Q estimators. 

We thank Niayesh Afshordi, Wayne Hu, Dragan Huterer, Bruce Winstein, and 
Matias Zaldarriaga for collaborative work and stimulating discussions. This 
work was supported by the Kavli Institure for Cosmological Physics through 
the grant NSF PHY-0114422. 
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